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I introduction 


The Luenberger observer Cl, 2, 3D allows extraction of the state of an 
observable linear system when given .1* the system input, '2 the system 
output, .3 the form of the system, and v4" the parameter values of the system. 

In those cases for which the system parameters are unknown the state observation 
is subject to error. Some previous i nvestlgators of parameter ignorance in 
observers C4,5D alleviate to some degree the observation error, but they are 
unable to guarantee the error vanishes or that their computational algorithm 
converges when the magnitude of parameter ignorance is large. We have previously 
reported C&3 the basics of a full order adaptive observer which negates these 
disadvantages. Our present paper considerably simplifies the exposition of 
the previous paper and extends, both computationally and theoretically, the topic 
of that paper. Briefly, the full order adaptive observer for single- input 
single-output observable continuous stable linear differential systems in the 
absence of a deterministic or random disturbance vector guarantees the vanishing 
of observation error regardless of the size of the constant or slowly varying 
parameter Ignorance. The observer parameters are directly changed in a 
Liapunov adaptive way so as to eventually yield the unknown full order Luenberger 
observer. The observer poles may throughout be placed freely in the stable 
region and no derivatives are required in the adaptive law. 

The Problem 

A differential system is assumed of the form 
w = Aw + Br w(0) = w° 

y = Cl 0 0 03w = Cw (I) 

A nxn 


B nxl 
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for whtch only the single output y = Cw = W| is available for measurement. 

It is assumed that a similarity transformation has been made if necessary so 
that the single- input *■ single-output system has fhe form of (I). It is 
assumed that seme or all of the elements of matrices A and B are unknown, A 
stable, w° may be unknown, and the pair (C,A) is completely observable. The 
observer is of the form 

z = Kz + GCw + Dr + Hu z(0) = z° 


i s 


( 2 ) 


F nxn G nxl 

D nxl H nxn and diagonal 

where K is arbitrary and u is a control vector yet to be defined but with the 

property that u->0 as t-*® . The problem is to adaptively form a triple (G,D,T) 

so that the error vector defined as e=z-T -l w vanishes as the system adapts. T 

is a non-singular square matrix with the property that CT=C. Fig. I illustrates 

A 

the adaptation. T(t) is a matrix which varies according to the adaptation 
procedure so that, when the adaptation is completed T(t) becomes T (i.e. 
lim T(t) = T) . 

t-H® 


Define a transformation x * T 'w so that e = z-x. Then (I) becomes 
x = A q x + T -, Br x(0) = T" 'w 0 

y = CTx = Cx 


( IA) 


and (2) becomes 


z = KqZ+ GCx + Dr + Hu 


z(0) = Z 


0 


(2A) 


It is desired for subsequent development that A n = T *AT be in the !: cutput 


form 
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a 1 1 I 0 0 . . . 0 

a£i 0 I 0 . . . 0 

| 0 0 i . . . 0 

• • • • t 

4 • • * • 

• • • • • 

a . 0 0 0 , . . 0 

nl 

wherefn the first column contains the system parameters and all other elements 
are zero save the super diagonal elements, which are unity. It is clear that for 
any non-zero matrix A with a single invarlent polynomial there corresponds a 
similar matrix Aq, although the elements of the similarity transf ormat ion may be 
unknown if elements of A are unknown. The following theorem defines the addition- 
al restriction that must be placed upon A so that both Aq = T *AT and CT = C. 
Theorem 

Cproof given in Appendix aD- Let A be an nxn matrix, C = [k,0,0, . . .Cf] a 

Q 

Ixn matrix with k^O, A Q an nxn matrix in output form, and T = C . • . 3 an 

T -| 

nxn nonsingluar matrix. There exists an (n-l)xn matrix T such that A=TA Q T iff 
the pair (C,A) is completely observable. 

As a result of the theorem, any observable system (I) may be placed by simil- 
arity transformation into system < I A ) with CT=C. The elements of T may be unknown 
since A is unknown. The problem will be considered as defined by equation ( I A) 
and (2A), so that e = z-x must vanish. Eventually the problem of constructing w 
from x wi I I be solved. 
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The Adaptive Law 

It Is now assumed, more for explanatory purpose than actual practical need, 
that some stable '‘nominal” plant matrix is either known ®r is chosen so that 
Aq = Aq + aAq, where Aq has all known elements and is in output form. Consequent- 
ly AAq contains all zero elements except for the left colurmwhich has elements 
that are to be adapted. Letting e = z-x, the vector error equation is 

e = Ke t (K+GC-A 0 -AA 0 )x + ABr + Hu 
e(0) = e° 

where AB = D-T *B. A theorem of Luenberger CG allows the eigenvalues of Aq~GC 
to be arbitrarily placed by selection of G (with the sole exception that Aq-GC 
cannot have the same eigenvalues as Aq). For the above error equation, let 
G=Gj t G^ and K = Aq-G^C. Then as a result of the theorem of Luenberger and 
of the special forms of A Q and C, the vector error equation is 

e = K Q e t (GjC-AAq)x + ABr + Hu (3) 

where K Q is an arbitrary stable constant matrix in output form with eigenvalues 
differing from Aq. The adaptive strategy is to change Gj and AB to eliminate 
the influence of x and r In (3); since by assumption Kq is a constant matrix, 
changing G| is equivalent to changing G and will be considered as such in the 
ensu i ng . 

For notational convenience in the next sections the following definitions 
are made. 
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H = diag CO, h h and nxn 

3 n-2' n-3 • 0 

n = order of plant 

(* 

m = number of zeroes in system transfer function. 

The error between plant state x and observer state z may be measured on!- 
by the scalar state variable G| = Zj-y = Z|-Xj. insure that only avail-.’ 
measurements are called for in the adaptive laws, (3) is "collapsed" to y *!d 
a scalar differential equation of the form 
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n 

I 

i=0 


I *A ]) = i l 


n- I n- 1 - j 

l l 

j=0 !=0 


l+J 


U) (j) 
a i+j x i 


m m-j 

+ I I 

j=0 1=0 


1+J 

I 


A 


e (1) JV 

e !+j r 


( 5 ) 


n-2 

+ l h j u 

1=0 


(1) 


I U 1 


For simplicity let Aj be a real characteristic value of A^-GC. Letting 
p=d/dt, the left side of (5) may be written as 


f n-l , 
(P + M I a. p 
1 *• i=0 


I 


where the a., 0 i <_ n-l, are defined by equating the above expression with 

the left part of (5). If It is desired to have no real observer pole, an 

obvious modification to Eq. 6 Is required. Now a reduction of order technique, 

similar to that of Gllbart and Monopoli [7U, Is applied to (5). The result is 
( n-l t n-l r n+m 


j” n-i j'j n-i j) f n+m 

(p+A.) I a . p e = l a.p l <j>.v 

'•i=o ' 1=0 ) L i=o 1 1 


n-2 

+ f + f + y h.u! 
x r j=0 J J 


(j) 


In which, assuming m < n-2, in (6) 


* f = S 


u 


. - a, a . 
i i n-i 

i = 0,1,2, 

n- 1 

i = n-l 

1-n 

i = n , n+ 1 , 


( 6 ) 


l .* 


and, defining the "state variable filters" v^. 
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n-l 

l a j V i 

j=0 


<J> .xi” 


= 0, l, 2, 


i-2 


( 8 ) 


= x, 


n-i 


i=0 


a . v 


<j> = r 


(i-n) 


I = n , n+l , n+2, • 


2(n-0 


= r 


v 2n- I 

"; 2 k ;' ai, * k 3 


f = I I ,.j 
T x k=l j=0 dt 


n _2 n-2 n-2-j 


J 


- i s L "w d+J 

k=0 J=0 I'O 


*,.3 


K * 


m iv- « r 

f > I l -J C *n +k 

r k -| j=0 d+ 


r (k-J-l>3 


(9) 


(10) 


n+m n-2 n-j 


- i-2 


sir Cv l J ' U 


-Xjio .So a ' +j+ ' ^ 

, th er, (7) should be changed to the extent tt> 

Should m=n-l tnen 

1 = n,n+l, n+2, • ••' 2n " 2 


♦l = 


S-,_ n - a l-n P n-I 


n-i 


] = 2n-i 


and UO) Is changed to 


n-2k-l d j rl r ”<-J- n 3 

f - I l 31 C *" +k 

r k= i j=0 dt 

2n-2 n-2 n-2-j _d^ £ v < 5 1 . 3 

- I I & 3 ' + j +l df> * 

k =n j=0 l u 


(7 .a) 
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It Is noted that (9) and (10) contain no other derivative of <j> k , 

0 i k < n+m, but According to the adaptive law (16), <f> is an available 

measurement. By careful manipulation of the state-space generation of (8), 
it is possible therefore to remove the totality of terms in (9) and (10) from 
(6). Generation of (8), it is pointed out, requires no derivatives of Xj. 

It is further noted that, since <j>^ is the change in parameters due to adapt- 
ation, as adaptation is completed 0 <k< n+m, and consequently I Im u.=0. 

t->oo ^ 

Table I gives the terms in u. for systems ranging from second to fourth 

J 

order when (8) has been chosen to be generated by a normal form. The table 
can be extended. . : ■>*. 


The implementation of u. as described reduces (6) to 

n- 1 . n- 1 . n+m 

(p+X )( l a. p )e . = ( l a. p 1 ) C I ♦, v] 

1 1=0 1 1 i =0 1 1=0 1 1 


(II) 


n- 1 


Taking Laplace transform of (II) and dividing by J a.s 1 yields 

1*0 ' 

a i .-o£( initial conditions) 

(s+A , )e ( = L ^ 4>, v f ] +^ 


1=0 


n-l . 

I a | s 

1=0 1 


( 12 ) 


for which fol lows 


n+m n- I 

®i + x i e i ~ I 4>i v i + I V, exp C-X.t] 
1 1 i=0 ' 1 1=2 ' 1 


( i 3) 


where are unknown constants or time dependent functions depending upon +:. 

n- 1 

initial conditions and {A.}, the set of characteristic values of j a. s'. 

1 1=0 1 
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A Liapunov function is now to be formed so that stability of the adaptive 
observer may be assured. To this end a positive definite function of the 
measured error e^ and the unknown parameter errors <j> j is defined as 


I 7 n+m 

V = (me; + l m *f) 

b 1 1=0 1 

Following Shackcloth [8], V can be made to be 

V = + e ( J ip } expC-X.t] 


(14) 


when 


m 


A = • S7 v i e i 


0 < j < n+m 


(15) 

(16) 


Other adaptive laws can easily be chosen instead if it is desired to increase 
convergence speed [9,103. 


Implementation of the adaptive law in (16) can be accomplished by 
enc 
example 


reference to (7) and to the definitions of the variables aj and For 


Vi ■ “n-i 


n- 1 


m 

s 

— x.e. 

Vi 1 1 


m 


V2 = “n-2 " 3 n-2 Vi = S n-2 + a n-2 m T . *1 e l 

n- 1 

m 

s 

v „e, 

m „ n-2 

n-2 


etc. 

in which may be ascertained. 

From the form of V, ej is stable in the sense of Lagrange with the re" ! or: 

of attraction determined by the unknown constants and the decaying exD'.'usntfa i 

time function. Clearly the region of attraction shrinks exponentially wvn 

time and eventually vanishes; consequently e ( is eventually asymptotical.'. 

stable and 1 1 m e . = 0. All derivatives of e. must vanish in the limit as xa i ' 
t-*» 1 1 
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since the scalar error equation (13) is linear and of first order and possessing 
finite frequencies. 

However, the Liapunov function (14) is defined on a non-compact manifold. 
Consequently {<}>.} is shown to be (eventually) stable but not necessarily 
asymptotically stable Cl ID- It is evident from (3) that each must vanish 
by adaptation In order to observe the correct plant state. Theorem 2 defines 
the restriction placed upon r(t) in order to guarantee vector error convergence. 

Theorem 2 CProof In Appendix Ef] 

Suppose there exists no set of real constants { q j > , i =0 ,1,2,..., n+m, for 
which the (observable) system command Input r(t) in its steady-state condition 
is a solution of the homogeneous differential equation 
n+m 

I q j r ‘ > =0 

i =0 1 

where n and m are defined in (4). 

Then I im <j>.(t) = 0, 1=0, 1,2, ..., n+m, and I im £(t) = 0 is assured. 

t-*» f-xc 

Coro 1 1 ary 

If the steady-state command input r(t) is periodic, a sufficient condition 

in order for lim _e(t) = 0 In (3) is that r(t) contain at least [n+m+l]/2 distinct 
t-*» 

frequencies in its steady-state condition. 

It is noted parenthetically that the corol I ary seems a generalization of a 
theorem of Lion Cl 2D although the applicability of that theorem to the present 
topic appears obscure. 

Reconstruction of T 

Using the ’‘nominal*' matrix Aq as Initial condition, the actual value or the 
system parameters may be determined by integrating the change In parameter:; 
defined in (16), until adaptation is complete and combining appropriately In the 



form of the matrix T. Thus T(t) ,: drifts i; toward T as adaptation progresses and 

'v 

I im T(t) = T. The example makes this technique clear. 
t-*“> 

A 

w, the estimate of w, is constructed from the observer output z by forming 

T(t)z. Consequently I im w = w. 

t->“ 

Practical Cons iderat ions 

Reference to Table i reveals that, speaking practically of analog implementa- 
tion, for high order systems a prohibitively large number of multipliers must be 

employed to generate the observer input u.(t). Since the magnitude of each u.(t) 

J J 

depends upon the magnitude of parameter change due solely to the adaptation 
process, it is reasonable to inquire whether u^.(t) can be omitted altogether 
(i.e. make H = 0) , espec ia I ly when the adaptation proceeds slowly by choice of const- 
ants in (16). By analysis of (14), it may be seen that omitting Uj(t) tends to 

degrade the adaptation process due to the inclusion of disturbances f and f , 

x r 

eq . (9,) and (10), 

However a theoretical analysis of a second order system Indicates that the 
u term may be safely omitted when the observer eigenvalues lie left of a curve 
passing through the left half-plane. This curve represents a trade off between 
frequency filtering in the adaptive law and magnitude of the adaptive gains. 
Generalization of this work awaits completion. 

Examp I e 

A third order plant with one zero Is considered for illustration. Let the 
plant be described by 



‘ 0 

i 

0 


' 0 

w = 

0 

0 

1 

w + 

C | 


; ( V V 

- ( a j +a | ) 

"‘W. 


. c o 


y = w. 
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In which cXq, a | Cq and Cj are unknown, a^, a^, a^ are the nominal 
I n output form, (I*) is 



-( a 2 +a 2 ^ 

1 

0 ' 



0 “ 

X = 

- ( a | ta | ) 

0 

1 

X 

+ 

b,-6| 


_ " (a 0 +a 0 ) 

0 

0 



\ 

0 
ca 

1 

O 

.Q 


y = x, = w, 

The error equation (3) is now 


CNI 

1 

1 

o' 


a 2 



0 



0 

- k l 

0 

1 

e + 

a l 

y 

+ 

e i 

r 

+ 

u l 

T 

1 

c? 

0 

0 


l 

P 

.? 



O 

1 



u o 

1- +* 


and the scalar error equation (5) Is now 

M ... . 4 
e l + k 2 e | + k |®| + k 0 e | = (a 0 + a | + a 2^ X | + ^ a l + ^ a 2^ x i + 2 X | 

+ B,r t (b q t 3 t )r t u ( t Uq 

Employing the definitions given In (7), (8), and Table I when n=3 and m : 
scalar error equation (5*) is equivalent to 

2 2 r 

(p+*|) Cp + a | p + a 0 )O| = (p + a ( p + a Q ) ( £ $ j v.) 

when Uj and u^ have been implemented as 

u 0 = + a l V 0 {,)) + 4> | ( v j ( 2 ) + a | v (D) 

+ d 3 (v 3 (2) + a j v 3 (l)) t d 4 (v 4 (2) + a ( V 4 (l)) 

“| “ i 0 V’ + *4 V | U> + *3 V 3 U) + *4 V 4 (l> 

As an illustration of the generation of Vj(j) appearing above, consider 


+ a , 


+ °0 v l 


" X, 


values . 

( IA*) 


(3*) 


(5*) 

I , the 

(II*) 


V| which 


is defined by 
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and the generation of v j ( j ) in normal form. 
Then, what is equivalent for v.. 


V ( 1 )' 


V t (2) 

ii 

_ 

L 


" a 0 _3i 



V (if 

+ 

1 ' 


v J (2) 


- 3 | 


.. _ 




Consequently Vj = v ^ ( I ) and both Vj(l) and V|(2) are available for measurement. 
Other v. (j) are generated in a similar manner. 


Defining cj>. as in (16), the observer has the form 


where 


z= 


-k, l 0" 


9 2 ‘ 


i 

o 


i 

0 

1 

i 

— o 

o o 

— o 

i i 

z + 

— o 

CD CD 

1 

y + 

b l 

b o 

"1 

+ 

1 

C C 

£ “ 


b. = - -2- 


b o = 


92 = 


9. * 


9 0 = 


m 

c 





e 

, v „ 


% 


1 4 


m 




s 





e 



m -. 


1 3 


3 




m 

c 




b 

e 

i x i 


rn 2 


1 1 




m 

m 

e , 


( — v . + a . 

— x. ) 

1 


m ( 1 1 

m^ 1 



m 

m 



s 

s 

e , 

( 

— v_. + a n 

— x i 

1 


m 0 0 0 

m 2 1 


and 


w = 


an - 


a, - 


9o d+ 


g, d+ 




3n - 


g 2 d + 


z=T(t)z 
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w is the estimate of plant state w, and I i m Note that CT = C. 

t-Xo 

A Simulation 

The third order system of the example was simulated on a digital computer 
using the following parameters 

a 0 


= 24 

o 

n 

o 

a 

C f = 30 

<3- 

(N 

II 

o 

m^/m^ = 

8000 

= 26 

a. = 74 

C 2 = 195 

k ( = 26 

m 0 /m 5 = 

2000 

= 9 

= 0 

o 

re 

II 

jQ 

k 2 = 9 

9 0 = 9 2 

= 0 


The eigenvalues of the observer (determined by {kj})were A ( = -4, A 2 = -2, A^ = -3. 
The input to the plant was a square wave of magnitude I and frequency 6t. Two 
parameters, b^ and g|, were adjusted by the adaptive law. These were initially 
set at b Q = 73, gj = -5 corresponding to a correct value of b Q = 75, gj = -74. 

Fig. 2 illustrates the behavior of b^, g^ e 2 , and e^ as a function of time. 

Remark 

A ^ A 

As has been previously stated, w = jz and I im w = w. In the general case of an 

* t-Xo 

arbitrary plant matrix A, the determinant of T may vanish for some instances of 

A 

time. These monentary occurrences of course, have no deterimental effect on w since 
* 

convergence of w to w is guaranteed, in the important particular case of the 
preceding example, however, advantage has been taken of the fact that det T is 

a | _ | 

constant by writing equation (*) as w = (T ) z. Since for the case of phase 
variable plant of high order the literal form of T * is easily produced, it is 

j _ | 

surmised that writing (T ) = T allows a particular simple construction of w 

when digital computation, rather than analog, is desired. 


Cone I us ion 

An adaptive observer has been demonstrated for single- Input single-output 
systems with constant or slowly varying parameters. Work is currently underway 
to extend the observer to multivariable systems as well as systems with rapidly 
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varying parameters and systems with noise. It is hoped that the adaptive observer 
will be eventually used not only for observing the state of an unknown system 
but in model reference paroblems and pole placement problems as well. 



- 16 - 


APPENDIX A 

In order to complete the proof of the theorem when A has repeated eigen- 
values (and the number of eigenvectors of A is less than n), the following 
lemma is needed. 

Lemma 

Let A be an nxn matrix, C = Ck, 0, 0, 0], k ^ 0, and \J»(A) = 

{all nonsingluar matrices Pjj = P ^AP) where J is an nxn matrix of Jordan form. 
Then there exists a Pe ip(A) such that CP = C C p j , 0^) , ( p j , 0_) D where 

each vector (p ,0) s [p . , 0, 0, ..., 0] has a dimension equal to the order of 
the corresponding i— Jordan block in J. 


Proof of theorem 


Since the similarity transformat ion matrix that transforms the matrix Aq 
( in output form) into the normal form A , 


A = 
n 


a . a _ ... a 

nl n2 nn 


is triangular, it suffices to show that the theorem is true with A Q replaced by 


(a) Let A n be partitioned into 


A' = 
n 


21 


M2 

\ 5 
72 
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where A j j scalar and is I x ( n- 1 ) row vector, Ajj and appropriately 
formed. Then 


A' = 


i 

0 


i < 

a I2 


k 

0 

J 1 

T 2 


_ A 2! 

A 22_ 



t; 1 . 


for I >_ I , and 

CA 1 = CTaJt' 1 = Ca't" 1 = C-A^ T^Tj * A^T" 1 ] 

since A|J = 0 for I <_ n-l. The term -Aji, T~*T| Is a scalar which may be any 
value. 


A 1 2 = C<5 |I 6 2l 6 31 ’•‘ <S (n-l)|- ] Where 


6 U = 


i=J 

‘0 I^j 


Let T = 


-2 


t 
— n- 


T 


where t. are row vectors. Then 
- 1 


Q, CC T ! A T C T: . A^C 1 !...; (CA n_ 1 ) ] = 


r k 

k. 

k 0 

’ k . 1 


1 

. 2 . 

. n- 1 

0 

* t, 

* t„ 

’ t . 

0 

• -1 

. -2 . 

. —n-l 

0 



• 

• 

• 
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Si nee nonsingular, t. mutually independent as well as independent of C. 
Therefore, observability matrix Q is of rank n. Therefore, A observable. 

(b) Suppose (A, C) completely observable but A J* T - *AT for any T. Let P 
be a modal matrix satisfying the lemma condition so that J = P ' AP and J is in 
Jordan canonical form. Let V be a Van der Monde matrix correspond i ng to J 
such. that A = VP 'aPV Since (A,C) observable, A possess but one invariant 
polynomial under supposition of the form of C [133; therefore V exists. By 
assumption^ f PV * . Since T is arbitrary, this implies that CT £ CPV or 
that 

CTV = CV ^ CP 

Because of the lemma, the P chosen in such that 

CP = [p,, 0 , 0 , ..., 03 . 

A necessary condition on CP for observability is that p ( be non-zero [143. 

By the form of the Van der Monde matrix, CV = [q|, 0, 0, ..., 03 with 
i 0 a constant. 

CV t CP implies q ( i Pj. But qj is any non-zero constant; thus, q ( r p ( 
implies p ( = 0 which Implies (A,C) not observable. Thus a contradiction is 


reached . 



Proof of Theorem 2 

It has been shown that 
I 1 m e (t) = 0 

t -x» * 
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APPENDIX B 


I im u.(t) = 0 for each j 

t-X» ^ 

I im <j>.(t) = constant for each ! 

t-H» 

Therefore from (7) and (7. a) 


1 im 

a. (t) = constant 

for 

each 1 

■f-xn 

1 



1 Im 

6 j (t) = constant 

for 

each 1 

t -^00 




Referring to n equations (3), each equation may be differentiated In 

so as to form the vector e = (t) = Cej n \ e^ n ..., e n 3 T . 

Employing ( B I ) and (B2) in determining I im e (t) and letting 3 . = 0 

+-w» “* • 1 


equations (B3) result. 


r, _ ( n- 1 ) 

0 - e„ + a 
2 n- 


y < n -l) + _ r <n-l> 


(n+l- i ) (n-i) . 

e. = e. , . + a . y 

I l+l n-i 1 


n- ! 
(n-i) 


+ 6 n _. f (n-I) I = 2,3,4, 


. .n- 


e n = V + 6 0 r 


All e’s may be easily eliminated from (B3) yielding 


( n- 1 


0 = 


l «, 




1=0 


m 


I s,* 

1=0 


Let the stable (but unknown) plant transfer function be 

c m 


f n 

l a j s 

M=0 


J 


y = 


l b l s 


i=0 


(Bl ) 


(B2) 


a manner 


for I > m. 


(B3) 


(B4) 


r 


(B5) 
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Comb ini ng (B4) and (B5) yields 

n- 1 . m . m . n , 

0 = C< l a : s ' ) (I b.s ) + ( l B.s 1 ) ( l a.s')> (B6) 

1=0 ' 1=0 ' 1=0 1 1=0 ! 

(B6) represents a condition upon r(t) which is assured in the limit, by.(gl) 
and (82) that is to say, after adaptation has forced e| to vanish. Two distinct 
possibilities exist regarding the solution of the (n+m)-th order linear 
homogeneous differential equation (B6): (a) either the steady-state system 
command input r(t) obeys (86) for some values ctj and or (b) the n+m+l co- 
efficients of polynomial in brackets are in the limit each zero. By supposi- 
tion of the theorem, (a) cannot occur; consequently (b) must be true. 

Using the assumption of cbservabi I i ty to insure that (B5) is relatively 
prime, it is easy to show by mathematical induction that condition (b) implies 
that the constants aj and B . are each zero, which was to be proved. 

Coro I I ary 

The corollary Is a direct result of placing the character i stic values 
of (B6) along the imaginary axis. The least number of distinct poles required 
of r so that it is not a solution of (B6) is exactly one more than the order 
of (B6), or n+m+l. Therefore r must contain at least n ^ — ■ distinct 
frequencies in its steady-state condition. 
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TABLE I 


Notation for generation of v. in normal form. In v.(j), j denotes 
the state variable, i denotes the function in (8). Example: for 

n=4, v | is def i ned by 


+ a 2 Vj 


+ a |V| 


+ 3 , 


v o = 


and generated by 


* 

V ( 1 )' 

- 

' 0 1 0 ' 


' v ( 1 )' 


0 

v, (2) 

= 

0 0 1 


v (2) 

+ 

1 

v. (3) 


_" a 0 " 3 I " a 2 


v ( (3) 


- a l 

- J 




- 


. 


n=2 

m=0 

n=2 
m= I 

n=3 

m=0 


n=3 

m=-| 


• * 

U 0 = ^2 V 2 + ^0 V 0 

U 0 <>0 V 0 + **2 V 2 

u 0 = •>! v | (2) + a | ^ J v | (l) + ^0 V 0 (2) + a l ^0 V 0 (1) 

• • 

+ a t $ 3 v 3 ( I ) + <j> 3 v^(2) 

= $ (Vj(2) + a j v j ( I ) ) + (v Q (2) + a, v Q < I ) ) 

+ (v 3 (2) + a | v 3 < I ) ) 

• • • 

u l = 0 0 v Q (l ) + pj Vj (I ) + b 3 v 3 (l ) 

u Q = (v q (2) + a, v 0 (l)> + (v j (2) + a, v ( (D) 

• « 

+ o 3 (v 3 (2) + a ( v 3 (l)) + <> 4 (v 4 (2) + a ( v 4 C I ) ) 



$ . 
i 


v.(l) 

i 


i^2 
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n=3 

m=2 


n=4 

m=0 


same as n=3 m=l 

with define in (7a) 

2 3 


u 0 ' J 0 ji, h a j V I (J) + *4 I, a j V 4<J> 


2 2 .2 
u i = ,J a j+i v i lJ) + ♦-> j, Vi v 4 (j) 


2 

U 2~ J 0 h v | (l) + ^4 V’ 


n=4 
m= I 


5 3 . 

u 0 ■ j 0 l 4 3 j v i < J> 

M3 


5 2 

u | = l l $j a i+ | v i ( J ) 

i =0 j = l 1 J ' 1 

M3 


u 2 = l $i v i (l) 
Z i =0 11 

i ?*3 


n=4 

m=2 


6 3 

u n = j j a . v. (j ) 

0 i =0 j=l 1 J 1 
M 3 


= 1 l <h 3 j+ | v j (j) 

1=0 j=l J 

i^3 


u = £ $ v (I) 

1=0 
M3 


n=4 same as n=4 m=2 

m=3 with (4> j } defined in (7a) 



PLANT 












